First-principles calculations of Cu(OOl) thin films: quantum size effect in surface 

energetics and surface chemical reactivities 
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First-principles calculations of Cu(OOl) free-standing thin films have been performed to investigate 
the oscillatory quantum size effects exhibited in surface energy, work function, atomic relaxation, 
and adsorption energy of the cesium adsorbate. The quantum well states have been shown and 
clarified at particular fc-points corresponding to the stationary extrema in bulk Brillouin zone, 
which are in good agreement with experimental observations. The calculated surface energetics and 
geometry relaxations are clearly featured by quantum oscillations as a function of the film thickness 
of the film with oscillation periods characterized by a superposition of long and short length scales. 
Furthermore, we have investigated Cs adsorption onto Cu(OOl) thin films as a function of the film 
thickness. Our systematic calculated results clearly show the large-amplitude quantum oscillations 
in adsorption energetics, which may be used to tailor catalysis, chemical reactions and other surface 
processes in nanostructured materials. 

PACS numbers: 73.61.-r, 73.20.At, 73.21.Ac, 



I. INTRODUCTION 



When the thickness of thin metal films approaches the 
nanoscale, the oscillatory quantum size effects (QSEs) 
associated with electronic confinement and interference 
will occur— i2i2i4 due to the splitting of the energy-level 
spectrum into subbands normal to the plane of the films, 
i.e., the quantum well (QW) states. These QW states 
lead to strongly modified physical properties and thus 
have been the subject of numerous experimental investi- 
gations in recent years 5 ^. For example, the QW states 
are found to be responsible for an unusual metallic film 
growth pattern 7 -^^, and for the thickness-dependent 
stability^! observed in the experiment. The QW states 
are directly connected to the oscillation in the exchange 
coupling between two magnetic materials across a non- 
magnetic spacer layer of various thickness 1 ^, and to giant 
magnetoresistanc o 13 ! 14 ' 15 ! 16 . Moreover, the QW states 
also give rise to an oscillatory phonon-electron coupling 
as the film thickness varies, and thus affect the transition 
temperature of the superconductivit y 17 i 18 . 

Experimentally, the characterization of the QW states 
are commonly measured using angle-resolved photoemis- 
sion spectroscopy (ARPES) and the scanning tunneling 
microscopy (STM). ARPES can be used to study the 
band structure along any direction of the surface Bril- 
louin zone (BZ), while STM offers the possibility to study 
local structures, such as islands, chains, dots, etc. Using 
the scanning tunneling spectroscopy (STS) technique, in 
which the differential conductance (dl/dV) is measured, 
one can reliably determine the energy of quantized elec- 
tronic states in the range of approximately 1 eV below 
and above the Fermi level. Theoretically, a number of ap- 
proaches have been used in the past in order to describe 
the electronic properties, in particular the QW states, 
in ultrathin metallic films. Quasi-one-dimensional mod- 
els, such as a square well potential— or the phase ac- 



cumulation model (PAM)2£, have been successfully used 
to interpret the energies of QW states. More sophisti- 
cated methods have also been used, such as the tight- 
binding approach 2 ^ and layer-Korringa-Kohn-Rostoker 
approac h 22 ' 23 . In a few systems the QW states have 
been investigated by self-consistent density-functional 
calculation o 24 i 25 i 26 i 27 i 28 i 29 i 30 ' 31 . There are strong reasons 
to use the ab initio methods. First, there are no ad- 
justable parameters, and a wide range of calculated struc- 
tural and electronic properties offer the possibility of a 
detailed comparison with experiments. Also, quantities 
such as the expected STM profiles and the amplitudes of 
the wave functions of the QW states, which cannot be 
obtained in simple approaches, can be calculated. 

In this paper we report our first-principles calcula- 
tions of the QSE in a specific QW system, i.e., the 
Cu(OOl) freestanding thin films. Previous QSE studies 
concerning Cu(OOl) are mainly focused on the oscilla- 
tions in magnetic interaction between the two ferromag- 
netic layers in fee M/Cu/fcc M (001) sandwich struc- 
tures where M denotes the ferromagnetic material. It has 
been demonstrated that the Cu(001) films have a long 
length scale of 5.6 monolayer (ML) and a short length 
scale of 2.6 ML of oscillation periods for magnetic cou- 
pling in the [001] direction, corresponding to spanning 
vectors at the belly and neck of the bulk Cu Fermi surface 
respectivel y 32 ' 33 ' 34 ' 35 . Recent experimental efforts have 
been focused on the new band structure properties of 
QW states in Cu(001) systen^^^^M 1 ^ 3 .^ 4 .. On 
the other side, the QSE of Cu(001) associated with its 
energetics was scarcely considered up to now. In partic- 
ular, there is no clear experimental or theoretic evidence 
of the interplay between the different oscillation periods 
of Cu(001) films by the belly and neck extrema in bulk Cu 
Fermi surface. In this paper we present a detailed first- 
principles study of the surface energetics of the Cu(001) 
free standing thin films. The QW states corresponding to 
the stationary extrema in bulk BZ are studied in detail. 
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The oscillations of the energetics versus the Cu(OOl) film 
thickness are identified and the corresponding oscillation 
periods are explained. We find that the quantum inter- 
ference of the QW states with different in-plane wave 
vectors result in a superposion between long- and short- 
length oscillating periods. 

The other purpose in this paper is to investigate the 
QSE character in surface adsorption energetics for a rep- 
resentative system in order to shed light on the effect 
of the QW states on the surface reactivities. Since the 
adsorption property is closely characterized by the chem- 
ical bonding between the adsorbate and the surface of the 
substrate, thus when the substrate is ultrathin, the QSE 
in the substrate will also influence the behavior of the 
surface adsorption. Here we address a particular adsor- 
bate system, i.e., Cs/Cu(001), as a case study to man- 
ifest the QSE in surface adsorption properties. Due to 
its simple electronic structure and active chemical prop- 
erties which is intrinsic for alkali metals, Cs is a unique 
adsorbate on metal surfaces, and has been extensively 
studied. Concerning Cs adsorption on Cu metal sur- 
faces, the recently published article a 45 i 46 ' 47 i 48 ' 49 i 50 mainly 
focus on the investigation of electronic, dynamic, and 
geometry properties of Cs layers on Cu(lll). The na- 
ture of interaction between the Cs adatom and Cu(OOl) 
film has been studied both via first-principles pseudopo- 
tential calculations^ 1 - or through and a jellium model 
approach^. 

Compared to the previous Cs/Cu(001) work, as men- 
tioned above, our present emphasis is on QSE in sur- 
face adsorption energetics, instead of giving a general 
overview of the nature of alkali-metal-atom adsorption 
on metal surfaces. In particular, although it has been 
well known that due to charge transfer from the adsorbed 
alkali atom to the substrate the alkali adatoms become 
partially charged^, the dependence of this charge trans- 
fer process on the thickness of the ultrathin substrate 
film is yet to be understood with high interest. We no- 
tice there are emerging some high-quality experimental 
data of QSE in charge transfer, via observation of the 
absorbed work function and photoemission as a function 
of the thickness of the ultrathin substrate film 54 . Fur- 
ther measurements related with the QSE effect of the 
adsorption are expected to be systematically reported af- 
terwards due to the availability of the high-quality quan- 
tum metal thin films. From this aspect a thorough the- 
oretical investigation is necessary and will be helpful for 
experimental references in near future. Our results show 
that in the ultrathin Cu(OOl) films, the surface adsorp- 
tion properties display well-defined QSEs. 

This paper is organized as follows: In Sec. II, the ab 
initio based method and computational details are out- 
lined. In Sec. Ill, the band structure and the proper- 
ties of QW states at the belly and neck points in bulk 
BZ are presented. In Sec. IV, the surface properties of 
the Cu(OOl) films, surface energy, work function, and in- 
terlayer relaxation, as a function of the thickness of the 
films, are presented and discussed. In Sec. V, the ad- 



sorption properties of alkali-metal Cs on Cu(OOl) surface 
is discussed in detail by presenting the sensitivity of the 
adsorption energy to the thickness of the Cu(OOl) films. 
Finally, Sec. IV contains a summary of the work and our 
conclusion. 



II. COMPUTATIONAL METHOD 

The calculations were carried out using the Vienna ab 
initio simulation package^ based on density-functional 
theory with PAW pseudopotentials— and plane waves. 
In the present film calculations, free-standing Cu(OOl) 
films by the so-called repeated slab geometries were em- 
ployed. This scheme consists in in the construction of a 
unit cell of an arbitrarily fixed number of atomic layers 
identical to that of the bulk in the plane of the surface 
(defining the two dimensional cell), but symmetrically 
terminated by an arbitrarily fixed number of empty lay- 
ers (the "vacuum") along the direction perpendicular to 
the surface. In the present study we have fixed the whole 
vacuum region equal to 20 A, which is found to be suffi- 
ciently convergent. The two dimensional unit cell of the 
fee Cu(100) surface is a square of edge with basis 

vectors a\ = o/2(l, 1, 0) and = a/2(l, 1, 0) where a is 
the Cu bulk lattice constant. The corresponding surface 
BZ is a square with two high-symmetry directions T-M 
and T-X. During our slab calculations the BZ integration 
was performed using the Monkhorse-Pack scheme^ 7 - with 
allxllxl /c-point grid, and the plane- wave energy cut- 
off was set 270 eV. Furthermore, the generalized gradient 
approximation (GGA) with PW-91 exchange-correlation 
potential has been employed with all atomic configura- 
tions fully relaxed. First the total energy of the bulk fee 
Cu was calculated to obtain the bulk lattice constants. 
The calculated lattice constant is a = 3.639 A, compa- 
rable with experimental^ values of 3.61 A, respectively. 
The use of larger fc-point meshes did not alter these values 
significantly. A Fermi broadening of 0.1 eV was chosen 
to smear the occupation of the bands around Ep by a 
finitc-T Fermi function and extrapolating to T — K. 



III. BAND STRUCTURE AND QUANTUM 
WELL STATES 

We first studied the properties of electronic struc- 
tures for bulk Cu and Cu(001) films. For this we have 
used two kinds unit cells of crystal lattice. Namely, 
unit cell I is defined via introducing the usual fee ba- 
sis vectors ai = (a/2)(j + k), a2 = (a/2)(i + k), and 
a 3 = ( a /2)(i + j), while unit cell II is defined by choos- 
ing Cu(001) as the basal plane, i.e., ai = (a/2)(i — j), 
a 2 = {a/2)(\ + j), and a 3 = ak. The volume of sec- 
ond unit cell is twice as that of the first one and it's 
convenient to extend to slab calculation since a 3 is nor- 
mal to Cu(001) surface. Correspondingly, the recipro- 
cal lattice basis vectors are bi = (2-k/o)(— i + j + k), 
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FIG. 1: (Color online) (a) GGA energy bands and density of 
electron states (right panel) of fee bulk Cu; (b) The Cu Fermi 
surface showing the belly and neck regions. The bulk and 
surface BZ are also depicted; (c) Calculated (GGA) energies 
at r in Cu(OOl) thin films as a function of thickness, with 
the energy set to zero at the Fermi level. The right panel 
plots the bulk energy dispersion in the [001] direction; (d) 
Calculated (GGA) energies at P in Cu(001) thin films as a 
function of thickness, with the energy set to zero at the Fermi 
level. The right panel plots the bulk energy dispersion in the 
[001] direction. The dotted line denotes Fermi level. 



b 2 = (27r/a)(i- j + k), and b 3 = (27r/o)(i+j - k) for 
unit cell I and t>i = (27r/a)(i — j), b 2 = (27r/a)(i + j), 
and b 3 = (27r/a)k for unit cell II. Figure 1(a) shows the 
band structure and the DOS of bulk fee Cu (unit cell I). 
The highly dispersive s-p band, typical of noble metals, 
is present. The more intense d band region is located 
between 2.0 and 4.5 eV below Fermi surface and can be 
measured by a large rise of the intensity in photoemis- 
sion experiments 12 . The s-d hybrid band region with 
two less intense dispersive states can also be seen in Fig. 
1 (a) . The Fermi surface of bulk Cu and the Brillouin zone 
for unit cell I are schematically shown in Fig. 1(b). Of 
particular interest for Cu(001) thin film studies are the 
two extremal points, i.e., the belly and neck points at 
the bulk Fermi surface. Previous extensive studies have 
shown that the stationary spanning vectors connecting 
the belly and neck points respectively play a key role in 
determining the oscillatory behavior of interlayer mag- 
netic exchange interaction Jm/m m fee M/Cu/M (001) 
sandwich structures where M denotes the ferromagnetic 
material. The belly extrema was found to correspond to 
long oscillating period in Jm/m while the neck extrama 
corresponds to short period. 

One inconvenience in discussing properties of Cu(001) 
film via the bulk Fermi surface calculated from unit cell 
I lies in the fact that none of three basis vectors of unit 
cell I is normal to (001) surface, which leads to an oblique 
projection of bulk BZ onto (001) surface BZ. On the con- 
trary, by use of unit cell II the discussions, particularly 
for QW states of (001) thin films, become very conve- 
nient. Beginning with this unit cell construction, the 
belly point in bulk BZ is still projected to V point, like 
with unit cell I. Whereas the neck point, after projection, 



turns to be at TX line. In this paper the neck point will 
be named P and P for bulk and surface BZ (with respect 
to unit cell II) respectively. Our band structure calcula- 
tion gives the position [in coordinates of (k x , k y , k z )\ of P 
to be (2tt/o)(0.4,0.4,0.0). 

Now we turn to study the electronic structure of 
Cu(001) thin film with focus on QW states. Since the 
previous studies 12 have shown that the long and short 
oscillating periods in interlayer magnetic exchange inter- 
action originate from band dispersions at T and P points 
respectively, we expect the other intrinsic film properties 
such as film energetics, atomic structure relaxations, and 
even surface chemistry are closely related with the QW 
states at these two kinds of stationary fc-points. Fig- 
ure 1(c) shows the energies of the QW states at T point 
as a function of the film thickness without interlayer re- 
laxation. The energy zero is set at the Fermi level of 
each film. Interlayer relaxation effect is also studied and 
it is found that the overall thickness dependence of the 
energies are similar to that without relaxation. For com- 
parison, also plotted in Fig. 1(c) (right panel) is the 
energy dispersion in the bulk along [001] (r — > X) direc- 
tion, which corresponds to the center T of (001) surface 
BZ and determines the energy range for the QW states 
shown in the left panel. One can see from Fig. 1(c) 
that as the thickness of the film increases, the energy 
of a given QW state also increases. When the thick- 
ness of the film is increased to be about 5.0 monolayers, 
then a quantum well state, with the energy crossing the 
Fermi level, occurs. The next energy crossing with the 
Fermi level occurs at the film thickness of ~11 mono- 
layers. Our calculated result of quantum well states is 
in good agreement with the recent experimental ARPES 
measurement 12 , in which the photoemission intensity as 
a function of energy for fixed Cu thickness and as a func- 
tion of Cu thickness for fixed energy clearly show the 
existence of the QW states of the sp electrons in the Cu 
layer. Fig. 1(d) shows (left panel) the energies of the QW 
states at P point as a function of film thickness without 
interlayer relaxation. The energy zero is set at the Fermi 
level of each film. Again the bulk energy dispersion along 
[001] direction (starting from P point). One can see that 
for every incremental increase in the film thickness of 
about 2.6 ML, then a new quantum well state crosses 
the Fermi level. A comparison between the left panel 
and right panel in Fig. 1(c) [or (d)] shows that the bulk 
electron band of Cu works very well for QW states as 
the Cu thickness is greater than 5 ML, which agrees well 
with the experimental results^. 

The quantitative understanding of the QW states 
showed in Figs. 1(c)- (d) is usually analyzed using the 
so-called phase accumulation model (PAM) 60 ' 61 . Here 
the free standing Cu(001) film is considered as a quan- 
tum well confining electrons between the two vacuums 
in the slab. Since the system is invariant on transla- 
tion parallel to the (001) plane, the in-plane wavevector 
k|| = (k x ,k y ) is a good quantum number. Thus, for a 
given k|| , the quantization condition for an electron state 
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in such a well is given by 



2k^Nd + 2(j> = 2im, 



(1) 



where N is the number of atomic layers in the film, 
d = a/2 the interlayer spacing, (j>g is the phase gain of 
the electron wave function upon reflection at the the film- 
vacuum interface, n is the number of half-wavelengths 
confined inside the QW and describes the Cu elec- 
tron wavevector component along [001] direction for the 
nth QW state. Equation (1) has been successfully used 
to model QW states in metal thin films, and its valid- 
ity has been rigorously tested by full-scale first-principles 
calculations for some special systems 2 ^,. It should be 
addressed that QW states will not only form at the T 
point but also form in a large part of the surface BZ, and 
the quantization condition (1) may be applied through- 
out the entire zone. Using Eq. (1) one can calculate 
the periodicity for the QW states crossing the Fermi 
level, AN = n/[kj,(\s.\\)d], where fcj;( k ||) is the bulk 
Fermi wavevector along [001] direction for a given in- 
plane wavevector k||. For the T point (kn = 0), from 
the right panel in Fig. 1(c), one can see that the upper 
branch of the bulk sp band runs through ~ 33% of the 
BZ, k{ = 0.33tt/c. One gets AiV = 6. Therefore a new 
QW state of k|| = (0,0) occurs every 5.6 ML increase of 
the film thickness, which is verified in the left panel in 
Fig. 1(c) that an energy branch moves down, crossing 
the Fermi level for every incremental increase in the film 
thickness of 6 layers. Similarly for P point correspond- 
ing to k|| = (27r/a)(0.4, 0.4) in bulk BZ, one can derive 

from the right panel in Fig. 1(c) fc^ = 0.757r/c. In this 
case, Eq. (1) gives the periodicity for the QW states to 
be AN = 2.6. Therefore, Figs. l(c)-(d) reveal that the 
periodic behavior of the Cu(001) QW states is essentially 
governed by the bulk Cu electronic structure properties. 
In particular, the extremal points of the Fermi surface 
and the extremal points of the QW branches close to the 



Fermi level occur at the same kn, i.e., at k|| = (0,0) in 
the bulk BZ corresponding to T in the (001) surface BZ 
and at the neck extrema k|| = (27r/a)(0.4, 0.4) both in 
the bulk BZ and in the surface BZ. This relationship be- 
tween the Fermi surface and the quantum- well dispersion 
turns out to be very fruitful for the understanding of the 
connection between the periods in energetics and surface 
reactivity discussed below. 

Using the energy data for QW states in Figs. l(c)-(d), 
and with the help of Eq. (1) and the bulk band struc- 
ture shown in the right panels in Figs. l(c)-(d), one can 
deduce the phase shifts <j> for the QW states at P and P 
points. They are plotted as filled circles in Fig. 2(a)-(b) 
respectively. Since the energy range of interest is quite 
far from the calculated vacuum level of 4.5 eV (relative to 
the Fermi level), this phase shift curve is nearly feature- 
less for most of the energy range, except for the feature 
near the top of the band. In general the phase shift de- 
pends on the incident (with respect to the film/vauum 
interface) energy and momentum of the electron. A sim- 
ple expression obtained from the WKB approximation is 
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FIG. 2: (Color online) The phase shift of the electronic wave 
function at the film- vacuum interface as a function of energy 
of QW states at (a) F and (b) P points. 



often used in the literature and has the form 62 
<P(E) = 



3.4 (eV) 



4.5 



(eV) - (E - h 2 kj/2m) 



(2) 



Here 4.5 eV is the Cu work function (in the thin thickness 
limit) , E is the electron energy measured from the Fermi 
level and kn is the the in-plane wavevector of QW. The 
image potential due to the electron-hole attraction is like 
a one-dimensional hydrogen atom. That is where the 3.4 
eV is from which actually equals one quarter of the hy- 
drogen ionization energy. For the belly of the Cu Fermi 
surface fcii = while for neck point we obtain fcii = 0.98 
A -1 . Equation (2) is plotted as a dashed curve in Fig. 2. 
Note that the data points within a given energy range 
can come from films of very different thickness. The 
datasets at both T and P vary smoothly as a function 
of energy with little scattering, indicating that Eq. (1) 
does determine the energies of the QW states quantita- 
tively, provided the phase shifts are known. 

The energies of QW states near the Fermi level can be 
measured by STS experiments, in which the local density 
of states (LDOS) is probed through the dl/dV curve. 
The distinctive sharp peaks in d//dU curve are charac- 
teristic of the QW states at different quantum numbers. 
The most easy-to-see quantity in STS measurement is the 
energy gap AE between the highest occupied QW state 
and lowest unoccupied QW state. To see the thickness 
dependence of AE, taking the derivative of Eq. (1) with 
respect to energy and evaluating it at the Fermi level for 
a given N, one has^ 
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FIG. 3: The inverse of the energy gap between the highest 
occupied QW state and the lowest unoccupied QW state as a 
function of Cu(OOl) film thickness. The squares in the figure 
are the results for F point in surface BZ, while the circles are 
for P point (projection of bulk Fermi surface's neck onto the 
(001) plane). 



where vf is the Fermi velocity obtained from the slope of 
the bulk band at the Fermi level and 4>'{Ef) the energy 
derivative of the interface electronic phase shift at the 
Fermi level. Therefore, the measured 1/AE curve is a 
linear function of N, with a slope connected to vp. Figure 
3 shows the calculated 1/AE at T and P points using 
the QW energies in Figs. l(c)-(d). One can see that 
the two curves follow a straight line with different slopes. 
Also it can be seen that whence a new branch of QW 
states crosses the Fermi level, then a kink occurs. Due to 
the different periods at T and P, the kinks in these two 
curves are also located at different values of N. Note that 
the intersection of the linear curves with the horizontal 
axis is not necessarily at N = 0, due to the nonzero 
energy derivative of the interfacial phase shifts at the 
Fermi level^i. 
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FIG. 4: Calculated electronic density of states at the Fermi 
energy as a function of Cu(001) thin film thickness. 
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FIG. 5: (a) Monolayer energy E(N)/N and (b) correspond- 
ing surface energy for fully relaxed Cu(001) lxl slabs as a 
function of thickness. The red curve is a least-squares fit to 
the surface energy. 



Another central quantity closely related with above 
discussed QW electronic structure of Cu(001) thin film 
is the film electronic DOS at the Fermi level, D(Ep). 
As shown in Fig. 4, the film D(Ep) with respect to the 
film thickness displays well-defined oscillations with os- 
cillation periods characterized by a superposition of long 
and short length. Note that in obtaining Fig. 4 we have 
increased the fc-grid for integration to 35 x 35 x 1 for 
each value of N, which ensures the precision of the re- 
sult. The oscillations in D(Ep) give an periodic band 
energy contribution to the total energy, thus producing 
oscillatory surface energetics and reactivity which will be 
shown below. 



IV. FILM ENERGETICS AND INTERLAYER 
RELAXATION 

Figure 3(a) shows the total energy per monolayer 
E(N)/N as a function of the number N of layers in the 
slab. The atoms in the slabs have been fully relaxed dur- 
ing the calculations. One can see from Fig. 3(a) that with 
increasing the thickness, E(N)/N gradually approaches 
a constant value which in the limit is equal to the energy 
per atom in bulk Cu. 

An energetic quantity more suitably tailored to QSE 
is the surface energy which is defined as one-half of the 
energy difference between the film and the bulk with the 
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FIG. 6: The second difference of the total energy as a function 
of thickness for freestanding Cu(OOl) thin films. 



same number of atoms, including the proper subtraction 
of a term linear in A^~. The thickness dependence of sur- 
face energy is calculated and shown in Fig. 3(b) (curve 
with squares) . It reveals that the surface energy approxi- 
mately follows a periodic oscillatory form, in accord with 
the oscillation properties of film D{Ep). This can be 
simply reasoned by the fact that the total band energy 
with respect to the Fermi level is related to the DOS by 
the equation E hand = J Bp (E - E F )D(E)dE. Therefore, 
each time the total DOS crosses the Fermi level, the band 
energy will also as a response have an oscillatory change, 
which consequently results in periodic oscillations in sur- 
face energy. In the same manner as discussed above, 
one can expect that the oscillations in the surface energy 
consist of a superposition of oscillations with the periods 
corresponding to the extremal points of the bulk Fermi 
surface. To illustrate this, also plotted in Fig. 3(b) (red 
curve) is a least-squares fit to the surface energy with the 
following expression 



i=l 



(4) 



where N is the number of monolayers (ML), Ai = 5.6 
and A2 = 2.7 are the periodicities (in unit of atomic ML) 
corresponding to the bell and neck extrema in bulk Cu 
Fermi surface, A4 and <pi are the fitting parameters. It 
can be seen that the fitting curve well reproduces the os- 
cillating behavior in the surface energy. Thus the periods 
extracted from the calculated Cu(OOl) surface energy are 
well consistent with the periods extracted from the bulk 
Fermi surface, and again the stationary extrema in Fermi 
surface play a key role in determining the oscillatory be- 
havior in surface energetics. 

Some authors have used the energies of QW states at 
L to discuss the stability of the film. See Ref^S for an 
example. The key point employed there is that even 
though the confinement takes place only in one of the 
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FIG. 7: (a) Work function for clean and Cs-adsorbed Cu(OOl) 
thin films as a function of thickness, (b) Planar-averaged 
electrostatic potential of clean (solid curve) and Cs-adsorbed 
(dotted curve) Cu(OOI) thin film (ten-layer Cu slab). 



three dimensions, the stability of the film could also be 
affected due to the variation in the electronic energy^. In 
the present case, since the above discussions have clearly 
shown that the neck point P in the bulk Fermi surface 
plays the same important role as the bell point (corre- 
sponding to r in surface BZ) does in determining the os- 
cillatory energetics properties of Cu(OOl) thin film, thus 
one can reasonably expect that the stability of the film 
also exhibits a superposed oscillatory behavior as a func- 
tion of the film thickness, with the oscillation periods 
consisting of long and short length. To illustrate this, 
following Refi^, we use the second difference of the total 
energy as a measure of the relative stability of an A^-layer 
film with respect to the films of N + 1 and N — 1 layers, 
which is defined as 

A(N) = E(N + 1) + E(N - 1) - 2E(N) (5) 

where E(N) is the calculated total energy of the fully re- 
laxed A^-layer film with the in-plane lattice constant fixed 
at the theoretical value. The result is shown in Fig. 5. 
A peak in the figure indicates a high relative stability for 
the film. It is not surprising that the stability is featured 
by a superposition of long- and short-period oscillations 
as is evident from the discussion above. Remarkably, one 
can see the N — 8 film is mostly stable while the N = 9 
film is particularly unstable. 

In addition to the surface energy, we have also calcu- 
lated the work function W of freestanding Cu(OOl) thin 
film. The work function, defined as the minimum en- 
ergy required to emit an electron from the surface to the 
vacuum, is one fundamental physical quantity for surface 
reactivity. An elementary picture of the work function in- 
volves a surface dipole layer that a valence electron must 
overcome in order to escape. Since the work function, as 
with many other properties, is also a function of the elec- 
tronic density, thus the changes in the electronic density 
by the crossings between the QW states and Fermi level 
will influence the work function in an oscillatory way. 
In addition to its relevance as an important element in 
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our understanding of surface science, a modified, or tun- 
able, work function can be useful for applications such as 
catalysis, because a slight change in the energy scale is 
exponentially amplified for chemical reactions^. Recent 
in situ experiments have measured atomic-layer-resolved 
work function and shown clear QSE in Ag/Fe(100) 66 
and Pb/Si(lll)££ systems. Here we have carried out 
first-principles calculations of the work function of the 
clean Cu(OOl) surface. The result is shown in Fig. 7(a) 
(filled squares) as a function of film thickness with re- 
laxed atomic geometry. One can see that the work func- 
tion is featured by an oscillatory behavior. As with the 
surface energy, the oscillations in the work function of 
Cu(OOl) thin films consist of a superposition of long- and 
short-length periods. 



TABLE I: Interlayer relaxations given in percent, Adi,i+i, of 
Cu(OOl) as a function of the thickness of the film. 



N 


Adi2 


Ad 23 


Ad 3 4 


Ad 4 s 


Ad 5 6 


Ad 67 


2 


-5.186 












3 


-3.351 


-3.350 










4 


-3.131 


+0.616 


-3.144 








5 


-3.425 


+0.081 


-0.099 


-3.431 






6 


-2.798 


+0.477 


-0.213 


+0.487 


-2.800 




7 


-4.076 


+0.126 


+0.398 


+0.394 


-0.115 


-4.066 


8 


-2.693 


+0.189 


-0.320 


-0.066 


-0.317 


+0.185 


9 


-2.917 


+0.369 


+0.133 


+0.484 +0.492 +0.134 


10 


-3.358 


+0.343 


-0.364 


+0.039 


-0.335 


-0.042 


11 


-2.805 


+0.415 


+0.256 


+0.415 


-0.029 


-0.025 


12 


-3.279 


+0.296 


-0.274 


+0.382 


-0.219 


-0.498 


13 


-3.699 


+0.237 


-0.322 


-0.134 


-0.407 


-0.938 


14 


-2.968 


+0.469 


+0.019 


+0.177 


-0.041 


+0.037 


15 


-3.286 


+0.127 


-0.315 


-0.245 


-0.721 


-1.015 



In addition to the above discussed film energetics, the 
relaxed atomic structures of the Cu(001) thin film will 
also be influenced by the occurrence of QW states at the 
Fermi level. Since the lxl supercells are employed and 
k|| is a good quantum number, thus only atomic inter- 
layer relaxation along [001] direction is allowed during 
our calculations. Here the interlayer relaxation, Adj^+i, 
is given in percent with respect to the unrclaxed inter- 
layer spacings, d , i.e., Ad hi+ i = 100(<i M+ i - d )/d . 
di^+i is the interlayer distance between two adjacent lay- 
ers parallel to the surface calculated by total energy min- 
imization, do = a/2 is the bulk interlayer distance along 
[001] direction. Obviously, the signs + and — of Ad^i+i 
indicate expansion and contraction of the interlayer spac- 
ings, respectively. The interlayer relaxations of Cu(001) 
films as a function of the film thickness is summarized in 
Table I. Furthermore, the interlayer relaxations are also 
plotted in Fig. 8 as a function of N for clear illustration. 
One can see: (i) The two outmost layers relax signifi- 
cantly from the bulk value, in agreement with the result 
from FLAPW calculation^. In the whole range of layers 
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FIG. 8: Interlayer relaxations of Cu(001) thin films as a func- 
tion of the film thickness. 



that we considered, the topmost interlayer relaxation is 
always inward (Adi^ < 0), with Adi,2 starting from the 
largest value of -5% for a slab with only two monolay- 
ers, and approaches a final value of -3% with increasing 
the thickness of Cu(001) films. Whereas, the second in- 
terlayer relaxation is always outward (Ag?2.3 > 0). Note 
that the first interlayer separation on most metal sur- 
faces is contracted, Cu(001) is one of the typical exam- 
ples, (ii) The interlayer spacings oscillate as a function of 
the thickness of the film with the period again consisting 
of long- and short-length scales. After 26 monolayers, 
which is the maximal layers considered here, the oscil- 
lations are invisible which suggest that the semi-infinite 
surface limit is now reached. 



V. ADSORPTION OF CESIUM: QSE IN 
SURFACE CESIATION 



In the above discussions we have extensively studied 
the QW states in Cu(001) thin films and the correspond- 
ing QSE in various physical quantities such as the surface 
energy, work function, and interlayer relaxation. To fur- 
ther illustrate the physical properties influenced by finite 
size of the thin films, in this section we focus our attention 
to the adsorption of Cs on Cu(001) thin films. Note that 
our study in this section is not just for an alternative ver- 
ification of QSE in solid thin film. On the contrast, the 
present study of surface adsorption as a function of the 
film thickness has its own cause to be emphasized. It is 
well known that metal surfaces are a prototype heteroge- 
neous catalyst, and have been widely studied in terms of 
the dissociative chemisorption of the reactant molecules. 
For a given solid, except for surface irregularities such as 
steps and defects, its surface reactivity is almost solely 




FIG. 9: The four different adsorption sites for Cs adatom on 
Cu(OOl) surface. 



determined by the crystallographic orientation. At re- 
duced size and/or dimensionality, particularly when the 
characterized size enters nanometer scale, however, the 
situation could be dramatically different from that of 
the bulk. In fact, size-dependent surface chemical activi- 
ties in metal films with the thickness of nanometer scale 
have been observed in previous experimental report o 69 ' 70 . 
Remarkably, in a very recent oxygen adsorption experi- 
ment on high-quality Pb(lll) thin films, Ma et al have 
clearly observed an oscillatory dependence of the chemi- 
cal reactivities on the film thickness^ 7 -, thus providing a 
most convincing proof of the key role played by the QW 
states in changing surface reactivity. It is expected that 
more experiments will address the issue of QSE in sur- 
face chemical reactivities. From this aspect, the present 
detailed theoretical analysis for the dependence of be- 
havior of Cs adsorption upon the Cu(OOl) film thickness 
is highly interesting and will be helpful for experimental 
reference in near future. 

Before we study the Cs adsorption properties as a func- 
tion of the thickness of Cu(OOl) thin films, we need to de- 
termine the energetically favorable adsorption site. Since 
the preference of adsorption site is not sensitive to the 
thickness of the substrate, thus to look for this pref- 
erence of the adsorption site, it is sufficient to give a 
study on the slabs with fixed thickness of the Cu(OOl) 
substrate, which at present is chosen to be 5 ML. We 
choose four most probable adsorption sites, which is enu- 
merated in Fig. 9. The binding energy is calculated 
using the following equation: Binding energy [atomic 
Cs]= -(£[Cs/Cu(001)j - £[Cu(001)] - 2£[Cs])/2 where 
£7[Cs/Cu(001)] is the total energy of a slab which includes 
two Cs atoms inside with a symmetric configuration, 
£'[Cu(001)] is total energy of the slab without Cs atoms, 
and £7[Cs] is total energy of a free Cs atom which is put in 
a 16 Ax 16 Ax 16 Asupercell. The calculated Cs/Cu(001) 
binding energies are 0.9189 eV, 0.9176 eV, 0.9174 eV, and 
0.9327 eV for site-1, site-2, site-3, and site-4 respectively. 
Thus site-4 is most stable for adsorption. This is in ac- 
cord with recent experiments that for substrate surfaces 
with square or rectangular symmetry the alkali atoms oc- 
cupy adsorption sites which maximize the coordination 
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FIG. 10: (Color online) (a) The contour plot of charge den- 
sity difference, p(Cs/Cu(001)) - p(Cu(001)) - p(Cs), in the 
(100) plane, for the cesiated Cu(001) surface. (b) (Left 
panel) band structure of cesiated Cu(001) thin film (5-layer 
Cu slab) with coverage O = 0.5; The right panel depicts the 
decomposition of local DOS for Cs adatom (dotted curves) 
and surface Cu atom (solid curves) into states with s (cyan 
curves), p (red curves), and d (black curves) character for the 
Cs/Cu(001)system. The inset in (b) pedicts the band charge 
densities around the Fermi level (lower panel) and above the 
Fermi energy (upper panel) of typical value of 3.5 eV. 

number to the substrate 7 ^ 7 ^ 7 ^ 7 ^ 7 ^! 7 !. Thus in the fol- 
lowing the atomic Cs is always put on site-4 during the 
simulation. This adsorption site is independent on the 
coverage O of Cs atoms. In this paper, we only consider 
the case of 9 = 0.5. 

To illustrate the surface bond formed between Cs 
adatom and Cu surface atoms, we depict in Fig. 10(a) 
the contour map of the electron-density difference for 
Cs/Cu(001). Here the plane depicted is normal to the 
surface and contains a Cs and two Cu's. The density dif- 
ference was obtained by subtracting the densities of non- 
interacting component systems, p(Cu(001))+p(Cs), from 
the density of the Cs/Cu(001) system, p(Cs/Cu(001)), 
while retaining the positions of the component systems 
at the same location as Cs/Cu(001). The solid and bro- 
ken lines denote an increase and decrease, respectively, 
in electron density upon Cs adsorption onto the surface. 
Covalent bonding is evident from the accumulation of the 
charge between the Cs adsorbate and the Cu(001) sub- 
strate. This charge is drawn principally from the atomic 
6s state of Cs adlayer. On the vacuum side of the over- 
layer the charge is depleted. Thus the Cs 6s state causes a 
polarization toward the Cu(001) surface. As a response, 
the work function will be decreased by this occurrence 
of surface polarization. On the other side, a significant 
counter-polarization of the Cs 5p semicore charge is also 
quite evident from Fig. 10(a). This effect opposes the 
reduction of the work function and the ultimate value 
of work function represents a self-consistent compromise 
between 6s and 5p polarization. The net result of these 
multiple surface dipoles, as shown in Fig. 7(b), is a low- 
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ering of the work function upon cesiation from 4.473 eV 
[clean ten-layer Cu(OOl) slab] to 2.259 eV for the cover- 
age 9 = 0.5, corresponding to the relaxed height of the 
Cs atom above the Cu film of 2.988 A. Also due to this 
opposite orientation between 6s and bp polarization, the 
work function of Cs/Cu(001) displays a non-monotonic 
variation with increasing the coverage of Cs: At low cov- 
erage, the work function decreases with increasing the 
coverage, going through a minimum (about at O = 0.5) 
and increases a little bit from there on to the high cover- 
age valued. The similar picture of cesiation process has 
largely discussed in Cs/W(001) system^,. A well-know 
conclusion in studying alkali-metal-atom chemisorption 
onto a metal surface is that although the region around 
the alkali-metal adatom is electrically neutral and no net 
charge transfer towards the metal surface is exhibited, 
the electrons of the alkali-metal adatom undergo a strong 
mixing with the substrate electron orbitals during the 
adsorption process. Such a process is local around the 
adatom and accompanied by a screening process, which 
is responsible for work function change. Furthermore, 
Figure 10(b) shows the band structure (left panel) and 
the orbital-resolved local DOS (right panel) for the Cs 
adatom and surface Cu atom, respectively. The complex 
Cs- induced charge rearrangement is more obvious. The 
Cs adatom experiences during adsorption a repulsive in- 
teraction of the valence electron with the induced image 
charge of the ionic nuclear core. As a result, the valence 
electrons shift upward in energy and hybridize with those 
of the substrate into bonding and antibonding state a 80 ' 81 . 
The bonding states are mainly around the Fermi level and 
are dominated by hybridization between adatom sp and 
surface Cu sp orbitals, as shown in the lower panel in 
the inset in Fig. 10(b). Whereas, the antibonding states 
shift up away from the Fermi level with a typical value of 
3.8 eV, and largely consists of hybridization between Cu 
s and Cs d orbitals, which can be seen from the upper 
panel in the inset in Fig. 10(b). A thorough descrip- 
tion of the bonding properties in Cs/Cu(001) system is 
beyond our intention in this paper. 

After finding the preferred atomic Cs adsorption site 
[site-4 in Fig. 9] and getting familiar with the Cs-Cu(OOl) 
bonding properties, we turn now to our central focus on 
the QSE in surface chemisorption and reactivity. For this 
purpose, we have given a series of calculations for the 
binding energy (the reverse of the adsorption energy) of 
the Cs adsorbate as a function of the thickness of Cu(001) 
thin films. Here as mentioned above we only consider 
the case of O = 0.5. The results are depicted in Fig. 
11(a). One can see that the binding energy of Cs onto 
Cu(001) thin films depends on the film thickness in a 
damped oscillatory way. These oscillations are featured 
by a superposition of long- and short-length periods, thus 
indicating a well-defined QSE in the surface chemical re- 
action of Cu(001) thin film with respect to Cs adsorption. 
In experiment this QSE of atomic adsorption can be ob- 
served by investigating the dependence of Cs coverage on 
the monolayers of Cu(001) thin films. We address here 
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FIG. 11: Calculated (a) binding energy of Cs adatom and (b) 
adsorbate height as a function thickness of Cu(001) films. 



that the calculated data in Fig. 11 have been carefully 
checked. This ensures that the oscillating behavior of 
the adsorption energy has nothing to do with the conver- 
gence problem, which may be encountered in the cluster 
calculations. The supercell approach retains 2D periodic 
boundary conditions during calculation, and can avoid 
this problem. Therefore, as with other quantities, the 
oscillations in adsorption energy is physically caused by 
the periodic change in the DOS at Ep. In fact, from 
what we have shown in Figs. 10, or from the simplest 
Anderson- Grimsley-Newns adsoptin mode l 82 ' 83 , one re- 
calls that as a Cs atom approaches to the Cu(001) sur- 
face (which can be approximated by a free-electron-likc 
metal), the valence s-p states of Cs adatom are broaden 
into resonant states due to as a result of interaction with 
the metal bands. Surface electrostatic potential will shift 
this hybridization down in energy to below Ep for sus- 
taining the whole neutralization. A lower DOS at Ep 
implies that the Cu(001) film has fewer electronic states 
to respond, whereas a higher D(Ep) means the film has 
more electronic states to respond, to the presence of the 
Cs adsorbate. Therefore, higher D(Ep) implies higher 
probability in the above hybridization process, causing 
the sp resonance to move to lower energies and to be 
occupied. This leads to a higher surface reactivity and 
higher adsorption energy. This physical picture can be 
formulated by expressing the adsorption energy as 

-Ead = -Eband(un-adsorbed) - E hlind (adsorbed) (6) 

r E F 

EADdE, 

where -Eband(un-adsorbed) is the band energy for a com- 
bined Cs/Cu(001) system with Cs adatom and Cu(001) 
film largely separated such that no chemsorption oc- 
curs, -Eband (adsorbed) is the band energy for the ad- 
sorbed Cs/Cu(001) system, AD =DOS(un-adsorbed)- 
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DOS(adsorbed) is the difference in the electronic DOS 
between the un-adsorbed and adsorbed Cs/Cu(001) sys- 
tems. As the film thickness changes, QW state en- 
ergy levels shift. Whenever a QW state crosses the 
Fermi level from above, it adds energy to the i^band (un- 
adsorbed) and adds electronic density to un-adsorbed 
D(Ep). Whereas Sband (adsorbed) decreases upon the 
crossing between the QW states and the Fermi level due 
to the above mentioned cause that the Cs-Cu orbital hy- 
bridization shifts more down in energy by the increase 
of D(E F ). Thus the net consequence is an anti-phase 
oscillation mode of the adsorption energy with respect 
the oscillation mode of D{Ep) as a function of the film 
thickness. 

The Cs adsorbate height is plotted in Fig. 7(b) as a 
function of the film thickness, which also shows the peri- 
odic oscillations indicative of QSE. Compared to the fea- 
tures in thickness dependence of binding energy, one can 
notice that oscillations in the adsorbate height arc much 
weak and the periods are somewhat not easily resolved. 
The cause may be due to the fully interlayer relaxations 
as a whole structure during calculations. If the substrate 
is not allowed to relax during the adsorption process, 
then it can be expected that there is a consistent corre- 
spondence in the oscillations between the binding energy 
and the height of the adsorbate. 

Moreover, we have also calculated the work function 
of Cs/Cu(001) as a function of the film thickness. The 
results are shown in Fig. 7 (open squares). Compared 
to the case of clean Cu(001) films, the work function is 
decreased by a typical value of 2.2 eV [see Fig. 7(b) for 
10-layer case]. Also, the adsorbed work function oscil- 
lates with increase of the film thickness. However, one 
can see that the oscillating mode is different from that 
without adsorption, and also the oscillating amplitude 
becomes much weak. This implies that the QW states 
are different for clean and cesiated Cu(001) thin films. 
As a result, there no longer exists a simple one-to-one 
correspondence between the adsorbed QW states and the 
bulk Cu Fermi surface. More detailed work will be done 
for this issue. 



dence of electronic structure, surface energetics, and in- 
terlayer relaxation upon the thickness of the film has been 
fully investigated, clearly showing the metallic QSEs of 
the film. These QSEs have been shown to be closely 
related with the occurrence of QW states at the Fermi 
level. Unlike some other simple sp metals such as Al, Mg, 
and Pb, the Fermi surface of bulk Cu is characterized by 
multiple stationary extrema such as the belly and the 
neck points. As a consequence, the Cu(001) QW states 
also display the abundant properties at these surface- 
projected fc-extrema. For example, we have shown the 
different film-thickness oscillation modes by calculating 
the energies of QW states in correspondence with these 
stationary extrema. Due to the interference between 
these two kinds QW states, the energetics and the cor- 
responding stability (Fig. 6) of the Cu(001) thin films 
have been shown in a consistent way to display a quan- 
tum beating behavior as a function of the film thickness, 
with the oscillation periods consisting of long and short 
length. In addition, the calculated energy gap between 
the highest occupied QW state and lowest unoccupied 
QW state, which can be directly measured via STS tech- 
nique, have been shown to display different slopes and 
kinks at these extrema in bulk BZ. 

We have also extensively discussed the other highly in- 
terested issue, i.e., the oscillatory QSE in surface chemi- 
cal reactivities, via studying the Cs adsorption onto the 
Cu(001) thin films. Through systematic calculations, we 
have shown that the Cu(001) surface chemical catalysis 
for Cs adsorption can be uniquely modulated and con- 
trolled by the occurrence of QW states at the Fermi level. 
The Cs adsorption energy and the work function after ce- 
siation are featured by quantum oscillations as a function 
of the film thickness. The amplitudes in these oscillations 
reach a typical value of 40 meV. From the application 
point of view, this result is quite intriguing because 1ML 
change in the film thickness can cause an effect equivalent 
to a temperature change as high as 400 degrees. These 
calculated results thus may be used as a guide to tailor 
catalysis, chemical reactions and other surface processes 
in nanostructured materials. 



VI. CONCLUSION 



In summary, the clean and cesiated Cu(001) thin films 
have been extensively studied by density-functional the- 
ory pseudopotential plane-wave calculations. The depen- 
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